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Abstract 



We develop and modify the Adomian decomposition method (ADecM) to work for a new type of non- 
linear matrix differential equations (MDE's) which arise in general relativity (GR) and possibly in other 
CJ , applications. The approach consists in modifying both the ADecM linear operator with highest order deriva- 

tive and ADecM polynomials. We specialize in the case of a 4x4 nonlinear MDE along with a scalar one 



5-H , describing stationary cylindrically symmetric metrics in quadratic 5-dimensional GR, derive some of their 

properties using ADecM and construct the most general unique power series solutions. However, because 
of the constraint imposed on the MDE by the scalar one, the series solutions terminate in closed forms 
exhausting aU possible solutions. 
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Q\ ■ 1 Introduction 
p 

Matrix differential equations (MDE's), which arise in many apphcations such as biomathematics, economics, 
^ ■ electronics, physics, stability analysis and statistics [1,2], are increasingly interesting fields of applied mathemat- 



5^ ■ ics and integral methods. In 5-dimensional general relativity (5DGR) and generally in J-dimensional general 
relativity (t/DGR) with d > A, linear and nonlinear 4x4 or 3x3 MDE's arise naturally in the 4 + {d — A) ox 
3 + {d — ^i) dimensional reductions, respectively [3-6]. Applications of MDE's to 4DGR are also known in the 
hterature [7, 8]. 

The kind of MDE's we deal with in GR are different from the Riccati, Sylvester or Lyapunov MDE's 
encountered in many other applications. In the latter MDE's, the coefficients of the unknown N x N matrix, 
denoted by are given matrices; while in the former case, these same coefficients are unknown scalar functions 
depending explicitly on the invariants of x and their derivatives. 

In order to introduce the MDE's encountered in GR we first define an abstract general form of them. Letting 
£i = ti{x'), 1 < i < N — I, and k = detx, we consider the following first order nonlinear MDE 

Gji£i,k,kpAp)ix'),p + Hj{£i,k,kpAp)x' =0, (1) 

where ", p" denotes the derivative d/dp with respect to the real variable p and Gj, Hj &Hq{1 < j <N —I) aie 
given explicit scalar functions of the unknown invariants of x and their first order derivatives. In equation ([T]), 
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X =1 denotes the N x N identity matrix and will be omitted throughout this paper. Notice that the powers of 
X higher than N — I and their traces do not appear in ([T]) for these can be reduced using the Cayley-Hamilton 
equation [9]. 

Introducing two arbitrary functions (Gq, T) and a matrix Q defined by 

Q^LGjX', (2) 

7=0 

we bring equation ([T]) to the form 

Q.p + TQ+Y.RjX^ = 0, (3) 

7=0 

where Rj = Hj - Gj,p - TGj for < j < - 1 . 

In many applications the extra functions {Go,T) depend, or are chosen to depend, linearly on ii{= trx)- 
Generally speaking, equation (O can be inverted to express as a polynomial in Q then substituted in ^ to 
bring the latter into a form where Q appears as the main variable. However, when the size of x is big such an 
inversion is not suitable; secondly, the method of resolution we are applying requires equation ^ to be kept in 
that form where X is the main variable. 

The aim of the paper is not to discuss in detail the properties of equation ([3]) or equation ([T|l, although some 
of them are formulated in section|2]and proven in sectionlH Rather, one of the purposes of this work is to develop 
and apply the Adomian decomposition method (ADecM) [10] to solve this type of nonlinear MDE's. ADecM 
is known as a reliable mathematical tool for solving algebraic, differential or integro-differential linear and 
nonlinear equations, but no general treatment for the type of nonlinear MDE's defined here is available. ADecM 
has undergone several modifications to extend its field of applicability, make it more efficient and speed up the 
convergence of the series solution [12-18]. In our application we specialize in the case of nonlinear MDE's 
describing stationary cylindrically symmetric (SCS) metrics (with four Killing vectors) in 5DGR extended by 
a Gauss-Bonnet (GB) term [3,4] and investigate them systematically. Knowing that no available technique of 
integration of differential equations can ever systematically handle this type of MDE's, our approach provides a 
systematic method for tackling them with the hope that it might be applied to MDE's in <iDGR. 

Some exact solutions to these MDE's describing SCS configurations in 5DGR have been constructed either 
by a perturbation approach [3] or analytically [4] and interpreted as cosmic strings. Since the MDE's are 
nonlinear, they may admit different solutions and the quest for further physically interesting exact solutions is 
still an open question, which may lead to construct cosmic strings endowed with similar geometric properties 
as the 4DGR ones. In fact, the 5DGR superconducting cosmic strings constructed so far [3,4] exhibit different 
properties than their homologs in 4DGR [19]. Besides, there has been recently a renewed interest in generalized 
a la GB GR theories [20-29] and particularly in string theory fueled by the discovery of Capodimonte-Sternberg- 
Lens Candidate, CSL-1, a pair of aligned galaxies which lie at a redshift of z = 0.46 whose double image could 
be a result of gravitational lensing caused by a cosmic string. 

At the time Ref [3] was written ADecM was being developed and not known to many workers (rather, it 
is still not known to many others), furthermore no other reliable method was available to tackle this type of 
MDE's, this is why we resorted to a perturbation approach. On the other hand, the alternative analytical method 
developed in [4] is not suitable for a systematic investigation of this type of MDE's. Thus, our aim is to further 
develop ADecM to work systematically for this type of MDE's in order to exhaust their set of possible solutions; 
this will turn out to be true for the case of the field equations describing SCS metrics in 5DGR extended by a 
GB term (sections [5] &[6ll. 

In section [2] we outline the properties of the MDE's ([T]) (or ^) and in section [3] we describe the physical 
problem and introduce the field equations — the Einstein-Gauss-Bonnet equations in the SCS case, which are of 
the form ([T]) with N = 4. In section |4] we develop and apply ADecM to the nonlinear field equations defined in 
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section [3] and prove some of their properties. The most general unique power series solutions are constructed 
in sections [5] & [6] then reduced to closed forms. These re-derived exact solutions are shown to exhaust the set 
of possible solutions to the field equations describing single SCS configurations in 5DGR extended by a GB 
term. Following the conclusion section, Appendix A completes the proof of corollary 3. and Appendix B that 
of Corollary 6. 

2 Properties of the matrix differential equations 

The formally defined nonlinear MDE (l3]l or ([T]) has the following properties 

1. Since the invariants of a matrix x remain invariant under a similarity transformation S, equation ^ 
remains invariant if X is replaced by 

X' = SXS'. (4) 

Since the rank of a matrix is also invariant under a similarity transformation [9], one can still generalize 
equation Q by making the functions {Gj,Rj), < j <N — \, depend exphcitly on the rank of X without 
affecting this property; 

2. equation ([3]) remains invariant under the simultaneous sign inversions X ~^ ~X & p ^ — p; 

3. It can be rigorously shown that any solution ;t(p) to equation (O commutes with its derivative x,p [30]: 

[X,X.p]=0; (5) 

4. Relying on ^ one can also show that any solution x{p) to equation ([3]) is a polynomial in a constant 
matrix A with scalar coefficients 

Xip) = a,ip)A' . (6) 

1=0 

Since Property 3. results from Property 4., the two properties are equivalent. It is also possible to first prove 
equation ^ using a decomposition method, as this is done in section |4] for the field equations describing line 
sources in 5DGR. 

3 The field equations describing 5-dimensional cosmic strings 

In 5DGR or Kaluza- Klein Theory (KKT) the spacetime has the topology of the product x S^, where is 
topologically equivalent to a 4-dimensional minkowskian spacetime and is a circle of radius a parameterized 
by the fifth coordinate x^: < < 2na. The fifth coordinate being spacelike and periodic, the spacetime has 
then one timelike coordinate and four spacelike coordinates. 

In 5DGR it is possible to modify the field equations by adding the GB term to the Lagrangian in such a 
way that the generalized Einstein-Gauss-Bonnet equations still involve at most second order derivatives of the 
metric [31]. The generalized action [32] reads upon ignoring a cosmological term, which is not relevant for the 
purpose of constructing cosmic string solutions 



3 



where 7 is a constant and the second term in the action is the quadratic GB contribution ££ = R Rabcd — 
AR^^Rab + R^ (the upper-case Roman indices take the values 1 to 5). The field equations derived upon varying 
the action with respect to the metric components gAs are then 



Rab - {l/2)RgAB + YLab = 0, (7) 

where Lab is the covariantly conserved {L^^ ■/, = ^) Lanczos tensor [32]: Lab = Ra^^^Rbcde — ^R^^ Raced — 
2 RacRb"^ + RRab - ( 1 /^)gAB ^. 

The 5-metric of the 5-spacetime reads 

d52 = g^5d/dx^, (A,B: 1^5). (8) 

We restrict ourselves to a SCS 5-metric which has four commuting Killing vectors {d(p,d^,dt,d5), df is 
timelike, and (5(^,^5) have closed orbits. Such a metric can be brought to the following form in an adapted 
coordinate system [3] 

d^^ = -dp2 + A«fo(p)dx"dx^ {a,b:2^5), (9) 

(the lower-case Roman indices take the values 2 to 5), where = p is a radial coordinate, = (p (periodic) & 
= z are the two other cylindrical coordinates, x^ = t (time), gn = — 1, gia = and gat = Kih{p) is a 4x4 real 

symmetrical matrix of signature ( 1- -). 

In GR the metric gAB, in our case A, is the main variable, however, in this paper, since we deal with the 
mathematics of the problem, we focus on the 4x4 real matrix = A^^A^p which we consider as the main 
variable. In terms of % the field equations of the pure KKT corresponding to the metric I© read [3,4] 

2Z,p + 4trx,p + (trz)Z + tr^^ + ip^xf = , (10) 
6trB + (tr;t:)2-tr;^2^0, (11) 

where B = x.p + (l/2)z^- Notice that (ITOl ) is of the form ([T) and ([TT]) is an extra constraint on the invariants of 
X- 

It is straightforward to bring equation (fTOl ) to the form ([3]). Subtracting the trace of equation (fTOl ) from 
equation (fTTI ) x 3 we obtain the equation tr;^^ = (tr;^)^ which we substitute into equation (flOl) to arrive at 

2(4tr;^ + 2x),p+tr;t;(4tr;t: + 2;t;) =0, withtr;^:^ = (tr;t;)^ (12) 

The system of equations (ITOl ) & (fTTI) . or its equivalent form (fT2l) . has been exactly solved in [3]. 
For a SCS 5-metric, the generalized Einstein-Gauss-Bonnet equations read [3,4] 

2Z,p + 4trXp + (trz)Z + ^X^ + {^xf + /{ (Z^),p " {^^X){X^),p 

+ [{^xf-^X^]X,p - (trz,p)te^ - (trz)z] - {l/2){iiX^),pX 

+ (l/2)[(tr;^);^3 _ (tr;^2)^2 _ {trx){Xrx')x + {^X?X\] =0, (13) 

6trB + {tixf- tr;t;2 + y[\r{Bx^) - tr (Bx)trZ 

+ (l/2)trB[(tr;^:)2-trz']}=0, (14) 

where the extra terms with respect to equations (fTOl) & ([TT]) are due to the Lanczos term. Equation (fT4l) is an 
extra constraint on the invariants of As we mentioned in the Introduction, some exact solutions to the system 
of equations ([T3] ) & (fT4)) have been constructed either by a perturbation approach or analytically [3,4] without, 
however, exhausting the set of possible solutions. As application of our developed ADecM we will investigate 
systematically the system of equations (fT3l) & ([14)) . Based mainly on ADecM and other techniques (differential 
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equations and matrix theory), our approach is cognitive leading to the unique exact solutions in closed forms, 
thus exhausting the set of all possible solutions to the system of equations ( fT3l ) & (fT4l) . 

The system (fT3] ) & (fT4l) is not readily solvable and need to be transformed to a reduced system as follows. 
Let the invariants of x be the functions /(p) = tr;^, g{p) = tix^, h{p) = trx^ and k{p) = detx, which have been 
called £i,^2)^3 and k in equation ([T]l- Let G{p) = g — and H{p) = h — f^. The Cayley-Hamilton equation for 
X reads then 

X' = fx' + {G/2)x^ + [{H/3) - {fG/2)]x -k. (15) 

We make use of the following relations which are valid whether the matrices X & X,p commute or not: 
tr(Z,pZ) = 5,p/2 & tr(;^p;^^) = /j,p/3. These relations together with the trace of ([T5] ) reduce ([141) to 

6[/p + {8/2)] - G + 7{(^p/3) + {fh/€) - ifg,p/2) 

- {fg/2) + (gV4) + (//12) - 2fc - (G/2) [/p + (g/2)] } = . (16) 

Now, the trace of ([T3] ). as provided by ([TT] ). also includes a term proportional to {2h^p + fh) as equation (fT6l ) 
does 

6/,p + 10/ + 8g + 7[2/j,p + //J - 3/g,p - %/p + 4//p - + / - = . (17) 

Hence, eliminating (2/j p + //j) by combining ([T6l ) and dlT] ). we obtain ^ in terms of f &. g: k = G(8 + 2/0 + 
3//^ +27/p). Consequently, any solution to the system (fT3] & [14)) is necessary a solution to the following 
reduced system where Q = 4f+{2 — yG)x — YfX^ + YX' 

2Q^p+fQ + 2G + jG^pX-7Gx^ = 0, (18) 
247^ = G(8 + 27G + 37/ + 27/p). (19) 

Equation dTSl ) has been readily obtained from ( [T3] ). In the case of pure KKT where the constant 7 = the 
system (HSl&illl) leads to 



4 The modified Adomian method 
4.1 Adomian 's method 

ADecM has been developed [10, 1 1] to deal with the initial and/or boundary value problems in nonlinear science, 
which can be modeled by 

LopM + /?op" + A'^opM = s , (20) 

where the unknown (scalar, vector or matrix) function u{x) is subject to some initial and/or boundary condi- 
tions. The function s{x) is a source term, (Lop +/?op) is a linear operator and Nop includes nonlinear operator 
terms. The method consists in splitting the linear operator (Lop + /?op) into two linear terms. Lop and R^p, where 
Lop represents the highest order derivative and is easily invertible and Rop groups the remaining lower order 
derivatives. 

Applying the inverse operator Lop' to both sides of ([201 ) leads to a Volterra integral equation 

u = w- Lop' (Ropu) - Lop' (A'^op") , 

where w is the sum of Lop' (s) and the terms arising from the application of the initial and/or boundary conditions 
to u. If one decomposes u into sum of components, u = Lm=o"m ' this leads to 

00 / 00 \ / 00 \ 

^ U,„ = W- Lop' ( /?op 52 ) - -^op' ( ^op ) • ^^^^ 

m=0 V m=0 / V m=0 / 
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Eq (I20I ) is solved upon fixing recursively the components m„,, which can be done by first spUtting the action of 
Nop on Lm=o"'n i'^ (I2TI ) into sum of terms ^^=0^™'^' where Af^'^ are called Adomian polynomials, in such a way 
that the first term Aq'^ depends only on uq, which will be the first component of u to be fixed, and A^'^ depends 
only on (mo, Ml), and A^'' depends only on (mq, "i , • • • , u,,,). Eq (|2T]) takes then the form 

00 / 00 \ / 00 \ 

£ M, = w - X /?opM„, - L-p' £ . (22) 

m=0 Vm=0 / Vm=0 / 

There is no unequivocal way to fix mq [17], however, a straightforward way to do it is to identify uo{x) with 
w{x) in (I22I) and the other components are determined by 

Um+i = -L^^{RopU,n)-L^^{A^), m>0, with uo = w, 

where A^** are evaluated using available formulas [10, 11, 18]. 

4.2 The modified method 

Without loss of generality, we consider equation (ITSl) and develop our method by extending ADecM to deal with 
this type of equations; the approach readily extends to the general case ©. We split x into sum of components 

00 

X=ZXn, (23) 

n=0 

resulting in 

CO 

/=£/„, with/,, = tr;t:„, (24) 

n=0 

and the linear operator Lop according to ADecM should be Lop = d/dp . However, it is sometimes fruitful and 
time-saving to adjust Lop to the kind of problem one is dealing with [12]. For the case of equation ([T8] ) our 
method consists to choose Lop = d/dp + (/o/2) the inverse operator of which reads 

L-pi[M] =exp(-Fo(p)) J\xpiFoip'))Mip')dp', (25) 
where M(p) is an arbitrary matrix and 

Fo{p) = {\/2)J' Mp')dp'. (26) 

Since we do not assume initial or boundary conditions, the action of L^p^ on the zero matrix, and only the zero 
matrix which appears in the r.h.s of (ITSl) . results in 

L-pi[0]=exp(-Fo(p))C, (27) 

where C is any constant matrix of integration. Replacing / by its expansion (l24l ) in equation dTS] ) then applying 
Lop* to the latter and using the fact that Lop^ [22,p + foQ] = 2Q, we obtain 

2Q = exp(-Lb)C - Lopi [(/-/o) 2 + 2G + yG.pX - yGx"] • (28) 



Next, we choose to split Q in the l.h.s of (1281 ) into a linear part Af + 2% and a nonlinear one Qnon = — jGx — 

Q = 4f + 2x + Qnon, (29) 
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and keep only the linear part on the l.h.s of (1281) . Now, we apply the standard approach of ADecM and replace 
all the nonlinear contributions in ( [28] ). that is 



!2„on and (£/„)e + 2G + 7G,pZ-7Gz' (30) 



by sums of components 



Qnon = XI 2«(Z0,---,Zn), (31) 

oo \ oo 

Y,fn)Q + 2G + YG,pX-YGx^=Y.B„{xo,...,Xn), (32) 

1=1 ' n=0 



in such a way that the components 2„ and B„ depend only on (xo^ ■■■,Xn) and do not depend on Xm for m> n. 
Our method consists in defining these components such that if A'^op denotes any nonlinear operator, say, one of 
the operators in ^ then Qo (resp. Bq) = Nop[Xo], Qi (resp. Bi) =Nop[Xo + Xi] -NopiXo], Qn (resp. B„) 

= ^op[ZoH ^ Xn] ~^op[ZoH ^Xn-i]- Hcncc, the components Qn and B„ are different from Adomian 

polynomials. For instance, the zeroth components are 

!2o(Zo) = -yGoXo - YfoXo + YXo > (33) 
Boixo) = 2Go + yGo,pXo - yG^xl , (34) 

where Go = tr^o - (i^Xo)^- With that said, substituting Eqs (l23]l, & ^ into (|28]l we bring it to 

oo oo / \ 

8/o + 4X0 + Z + 4;^") = exp(-Fo)C - £ 2!2„ + L^^' [B„] , 

n=l n=0 V / 

where the components 8/„ +4Xn, for « > 0, are determined recursively bjQ 

8/() + 4zo = exp(-Fo)C, (35) 
8/„+i + 4Xn+i = -2Qn - L-pi [Bn] , « > . (36) 

In order to derive an expression for Xn+i' for « > 0, we first trace ( [36l ) and use the fact that tr/„+i = 4/„_|_i 
(recall we omit to write the identity matrix /) to obtain the formula 

= -{2tra, + tr(L-pi[B„])}/36, (37) 



which we substitute back into (1361 ) leading to 

Xn+i = ^ + . (38) 

In the following two sections we shall discuss separately the cases / / & / = 0. 



'There is no unequivocal way to fix the zeroth component [17] nor an unequivocal way to fix the other components. So in OSb . we 
could also choose 8/o +4xo = exp(— Fo)C — 2Q() without, however, affecting the following discussion. 
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5 Nontraceless solutions: / ^ 

First we focus on the case / 7^ 0. Without loss of generality, we can always assume /o / 0; if this were not the 
case we could permute the indices of the components /„ and those of Xn so that to make /o / 0. The component 
Xo is derived in the same way as we did with Xn+i'- tracing ([35]) we obtain 



36/o = exp(-Fo)trC, (39) 
with trC / since /o / 0. We set C = 8(A + 2) then differentiate (HH) using (Eg]) to obtain 

trC 

36/o,p = -(/o/2)exp(-Fo)8(ttA + 8) 
= -(/o/2)x36/o 

leading to /o.p +/q /2 = the solution of which, introducing a constant of integration po, reads 

/o = 2/(p-po). (40) 

Substituting (gOjl into (US, we obtain 

exp(Fo) = p-po. (41) 
With /o and exp(Fo) given by ( |40l ) & (|4TI ). equation ( [39l ) reduces to trA = 1 and equation ( [35l ) to 

;t;o = 2A/(p - po) , with trA = 1 . (42) 

Knowing Xq & /o we derive the formula of the operator L^^ using equation (1251 ) 

^op' M = ^7-^ l\p'- Po)M{p')dp' . (43) 
P — Po J 

Next we determine the first component Xi- With Xo given by (l42l ) we obtain Go = 4(trA^ — l)/(p — po)^ 
which leads, using & to 



(A3-A2)-(trA2-l)A 

ip-poY 

8(trA2-l) (trA2-l)A + (trA2-l)A2 

^" = 7^^"^'" • 

With Lop' given by (|43] ). we obtain using (1381 ) 

47[2 (tr A^ - tr A^ ) + (tr A^ - 1 )2] 2 (tr A^ - 1 ) In | p - po | 



9(p-Po)3 9(p-po) 
4y{A^ -A^) + 27(tr A^ - l ) (a^ - A) 

(p-pay 



(44) 



Notice that & Xi already of the form ^ and all the other components are of the same form Thus 
on its domain of convergence , the matrix x = 117=oXn too is of the form In fact, if we try to obtain Xi we 
first derive Qi & Bi using the formulas 

Qi{XQ,Xi) = -riGo + Gi){xo + Xi)- 7(/o + fi){Xo + Xi? + YiXo + Xi?-Qo{Xo), (45) 
B,{xo,Xi)=2{Go + G,) + r{Go + G,).p{xo + Xi)- 7(Go + d ) (Zo + )' - Bo (Zo) , (46) 
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where Gi = tr {xo + Zi)'^ ~ i^^Xo + ^^Zi)^ — Go and Qq & Bq have been evaluated earher. Since both matrices 
Xo & Xi depend on the same matrix A (they commute), the new matrices Qi & Bi too are functions of the same 
matrix A, consequently, Xi given by (1381 ) in terms of Q\ & Bi is also a function of A which, after eliminating the 
powers of A higher than 3 using the Cayley-Hamiltoitl equation for A, reduces to the form i^. This argument 
applies to all remaining components. The whole argument remains valid and applies to the case / = 0, too, as 
we will see in the next section. 

Corollary 1. 

On its domain of convergence, the series representation, L^=o/t«' ^ solution to ( fTSl ) converges to a function 
;f (p) of the form Q. 

Notice that the constraint ( [T9l ) has not been used yet. Our procedure consists first in determining the com- 
ponents Xn, as we have done in (|42)) & (|44l) . and once these are known we evaluate separately the invariants of 
X (k, f & G) and substitute them into ( fT9l ). which turns into an algebraic equation. From this point of view 
equation ([T9l ) splits as 

247det(;^:o + Zi + ---) = 8(Go + Gi + ...) + 27(Go + Gi + ...)2 

+ 37(Go + Gi + ...)(/o + /i + ...)2 

+ 27(Go + Gi + ...)(/o + /i + ...),p. (47) 

Using (l42l ) & (l44l ) it is straightforward to check that the r.h.s of (|47] ) includes a term proportional to l/(p — 
poY emanating from Gq = 4(trA^ — l)/(p — po)^ while the l.h.s. does not. Consequently, in order to satisfy 
equation ( [T9l ) we have to choose 

Go = ^ trA^ = 1 . 
This new constraint on trA^ reduces the expression of Xi, equation (l44l) . to 

_ 87(trA3-l) Ay{A^-A^) 

9(p-po)3 (p-po)3 ' ^^'^ 

where the terms proportional to ln(p — po) no longer appear and consequently will not appear in the remain- 
ing components {Xnn ^ ^ 2). Hence, each component Xn^ n > 0, is a polynomial in x = l/(p — po) with 
matrix coefficients, which are functions of the same constant matrix A, thus they are commuting matrices. 
Since the two nonlinear terms in ( [30l ) remain invariant under the simultaneous sign inversions x & 
(p — Po) ^ — (p — Po) (Property 2.) and since both X'^ and X\ involve only odd powers of x, by (1381 ) all the 
remaining components ;^,„, ?w > 2, will involve only odd powers with the lowest power being 3; ;^o is the only 
component involving a factor of x. If we group the same powers of x^'+' , the general expression of X solution to 
the system dTS] ) & ( fT9b is a power series in x. This rearrangement or grouping of the terms of the series Y^^=q Xn 
is justified as follows. Each element Xab (a, ^ : 2 — > 5) of the matrix series is a power series in x which converges 
absolutely within some radius: \x\ < rah- If ''min denotes the minimum of all the radii rah, the series YJ^^^qXu 
converges absolutely in the domain \x\ < r^nm where the rearrangement of the terms is justified. 

Corollary 2. 

Any nontraceless solution to the system ([TSl l & ( fT9l ) has two series representations, a) Adomian series Y.7=oXn 
and b) a power series in x^'+' with constant commuting matrix coefficients M2i+i, which are functions of the 

^The powers of A higher than 3 result from the terms such as XoXu xj' stc in J45b & i46\ and are eliminated using dlSt where we 
have to replace x and its invariants by A and its invariants. 
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same constant matrix A subject to trA = trA^ = 1, of the form 



X = 2Ax + £m2,-+i (A);c2'+i , trA = trA^ = 1 , |x| < r^^ , (49) 

!=1 

where r^in is the radius of absolute convergence of both series representations. 

By Corollary 2., the other invariants of A, that is, trA^, detA & rank (A) are still free parameters. Corollary 
3. puts more constraints on these remaining invariants. 

Corollary 3. 

Any nontraceless solution to the system ([TSl l & ([T9l ) has a power series representation (|49l ) where the constant 
matrix A is subject to 

trA = trA^ = trA^ = 1 & detA = . (50) 

A proof of Corollary 3. is detailed in Appendix A. 

Now we proceed to fully determine the remaining components Xn (« > 2) using the recursive formula ( [36l ). 
Taking into account the constraints (l50l ). the expressions for Xo & Xi (equations (l42l ) & (|48] )) simphfy greatly 
and the Adomian series reduces to 

oo 

X = 2Ax- 47(A^ -A^)x^ + Y.Xn,\x\< r^,„ . (51) 

n=2 

The constraints dSOl l reduce the Cay ley-Hamilton equation for A to A^ = A^. Using this last equation in ( |45l ) & 
(l46l ) leads to 2i = & Bi =0, which we substitute into ([361) read^ 

8/2 +4;t2 = ''=^^2 = ^ Z2 = 0. 

With Xi = all the corresponding components are zero: G2 = 0, ^2 = & B2 = leading again to /3 = & 
X3 =0 and so on = for n > 2. Equation (|5TI ) reduces to the truncated solution with 

'"mill — °° 

P-Po (P-Po)^ 

where A is subject to ( l50l ). It is straightforward to check that ( l52l ) solves the system dTS] ) & ( [T9b . 

Corollary 4. 

Any nontraceless solution to (fTSl) & ( [T9l ) is a singular matrix (l52l ) with / = 2/(p — po), ^ = 0, trA = trA^ = 



trA^ = 1, detA = 0& r, 



mm 



The solution (1521) was derived earlier by a perturbation approach [3] then analytically [4] , but it has not been 
acknowledged in the literature that it is the unique nontraceless solution to (fTSl ) & (fT9l ). Other locally equivalent 
as well as globally nonequivalent solutions to (l52l ) can be derived by similarity transformations (HJl [3,4]. A 
classification of the solutions (l52l) according to rank (A) has been discussed in [4,30]. 

In the standard Kaluza- Klein 4 + 1 dimensional reduction, the components (§15/^55 = 0,^25/^55, A35/A55, 
^45/'^55) of the metric Q are proportional the electromagnetic potentials (Ai,A2,A3,A4), respectively, with A4 
being proportional to the electric potential. Rewrite Eq ([7]l as 

R^B - {l/2)R5^B = StiGT^sb , (53) 



We take L™,' [Bi] = L™,' [0] = 0: only one global constant of integration C is introduced in J27t . 
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where the r.h.s is the Lanczos term — {y/^hG)L^b reinterpreted as a 5-dimensional effective energy-momentum 
tensor source (here, G is the gravitational constant). For the solution (|52l ) the evaluation the I.h.s of (|53] ) depends 
exphcitly on rank (A). For rank (A) = 1 or 2, the I.h.s of (1531 ) includes only distributional contributions while 
for rank (A) = 3 it includes distributional contributions as well as continuous ones, which are proportional to 
7 [3,4]. The solution (l52l ) with rank (A) = 1 or 2 has been interpreted as neutral or charged cosmic string and 
that with rank (A) = 3 as an extended superconducting cosmic string (the continuous contribution of rJj-5 oo 
as p ^0) surrounding a naked electrically charged cosmic string core in longitudinal translation carrying [4] or 
not carrying [3] an electric current. 



6 Traceless solutions: / = 

Next, we consider the other alternative: / = 0. Notice that the above general formula, equation (|49l ). only 
applies to the case / 7^ since the trace of its r.h.s never vanishes identically for any values of trM2,+i (/ > 1). 
To analyze the case / = by the decomposition method we have to assume that at least one of the traces /„ is 
zero, which we take it to be /q. Hence, with /o = 0, equation (|26l ) leads to Fq = constant and the operator L^^, 
equation (1251 ). reduces to 

L-4[M] = fM[p')dp'. 

First we se^^ exp(— Fo)C = 4A in equation ( [35] ). with trA = since /q = 0, then apply the steps from ([35] ) to ([38] ) 
to obtain Xn (n > 0) as polynomials in p with commuting matrix coefficients, which are functions of the same 
matrix A. The first two components rather written in the form (O read 

_ _ 27trA^ - QtrA^ + trA^{S - YtrA^)p YtrA^{2 + p)A - lyA^ 

Zo-A, Zi- + 4 • 

It is obvious that x will have the form Within its radius of absolute convergence, i^mint we can rearrange the 
terms of X = lL7=oXn by grouping the same powers of p, the matrix x takes the following shape where trM, = 
(/ > 0) to ensure that / = tr = and where M/(A) commute since they are functions of the same matrix A 

00 

X = J^Mip' for p < rmin , with trM,- = for / > . (54) 

i=0 

Corollary 5. 

Any traceless solution to the system ([TSt & ([T9] ) has two series representations, a) Adomian series l^~=o and 
b) a power series in p' with constant commuting matrix coefficients M,, which are functions of the same constant 
matrix A subject to trA = 0, of the form ( [54] ). 

Substituting the series expansion ([54] ) into equations ([TSt . ( [T9] ) & ([T5] ) we bring them to the following 
expansions 

00 00 00 

(Ell)^ £c„p" = 0, (dSl)^ £c„p" = 0, (E])^ ££„p" = 0, (55) 

n=0 11=0 n=l 

where C„ & E„ are matrix coefficients and c„ are scalar ones depending on the traceless and commuting matrices 
M,. Setting these matrix and scalar coefficients to zero will alow us to determine the matrices M, (trM, = for 
/ > 0). 

In the calculation below we will make use of the Cayley-Hamilton equation for Mq (trMo = 0) 

= (trM^ /2)M(^ + (trM(] /3)Mq - detMo . (56) 



^The matrix A introduced for this case / = has nothing to do with previous notations. 
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The equations Co = 0, cq = & ci = lead respectively to 

(4 - 27trM^ + 6yMI)Mi + 1\xMI - /(tr M^)M^ 



-27tr(MoMi)Mo = 0, (57) 
detMo = trMl/{3y) + (trM2)Vl2 , (58) 
TtrM^tr (MqMi ) - 4tr (MqMi ) - 67tr (M^Mi ) = . (59) 



Adding (159]| to the trace of (l57])xM() leads to 

trM^ [trM^ + 3tr (MqMi )] = . (60) 



First using the Cayley-Hamilton equation for Mq, equation (1561 ). to eliminate any power of Mq higher than 3 in 
the product (|57])xMq then combining the trace of (|57] ) with that of (|57])xMq to eliminate tr (MqMi) leads to 

trM^(trM^-2/7) =0. 

From this last equation we have either 1) trMg = 2// or 2) trMg = 0. 



6.1 trM^ = 2/r 

Equation (1581 ) implies detMo = 1/7^- Using the Cayley-Hamilton equation for Mo, equation (l56l) . to eliminate 
any power of Mq higher than 3 in the product (|57])xMq. then combining the trace of ([571)xMq with ( [59l ) to 
eliminate tr (MqMi) we obtain 5trMg +6tr (MoMi) = 0. Solving this last equation together with equation (l60l) . 
which now reads trM^ + 3tr (MqMi ) = 0, results in trM,^ = & tr (MoMi ) = 0. Substituting these values into ( |57] ) 
reads 67MqMi + (4//) — 2Mq = 0, then multiplying it by /(I — /Mq), which is the inverse matrix of Mq, leads 
to 

Ml = -1/(37) + 2MoV3. (61) 

With this value of M\ , equation £1 = is identically satisfied. 

Equations Ci = & C2 = are brought to the following forms, respectively, using (|6TI ) 

[2(7 - 1 ) + 3y^\.v (MoM2)]Mo + if Ml - 9y^mIM2 = , (62) 
37^tr(MoM2)+97^tr(M(^M2)-2-97 = 0. (63) 



Solving (|63]) along with the trace of (l62l)xMn we arrive at 

tr(MoM2) = 1/7^ + 3/(272), tr(Mo'M2) = -l/(974) + l/(273). (64) 

Using the Cayley-Hamilton equation for Mo, equation (l56l ). to eliminate any power of Mo higher than 3 in the 
product (l62l ) xMq, then ehminating tr (M0M2) & tr (MqM2) by (l64l ) we obtain 7 = 2/3. With this value of 7 and 
the relations (l64l ). the matrix equation £"2 = reads 

(64/9)M^M2 - (16/3)MoM2 - (32/3)M(^ + 14 = 0. (65) 

Tracing both sides of (1651 ) and using (l64l ) along with trMg = 2/7 (7 = 2/3), we obtain —4 = 0. Hence, the case 
trMp = 2/7 does not lead to traceless solutions to the system (ITSl) & ( [T9l) . 
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6.2 trMj = 

Equation (1581 ) implies detMo = 0. In this case equation (l60l ) is satisfied so that trMg remains undetermined. 
Now, if we combine the ti^ace of (|57])xMq with ( [59l ) we obtain tr(MoMi) = 0, which reduces equation ( [59l ) to 
tr(M,)Mi) = 0. Substituting trMg = & tr (MqMi) = into ([57]), the remaining equation (4 + 67M,))Mi = 
leads upon multiplying by (1/4) — S/Mq/S, which is the inverse matrix of 4 + 67Mq, to 

Ml =0. 

Hence, the traceless case (tr Mi = for / > 0) with trM^ = leads to detMo = 0, Mj = and trM^ & Mq 
remain undetermined. As will be shown in Appendix B, even though trMp is still a free parameter, all the 
remaining matrices M,, / > 2, are proven to be zero. We have thus obtained the following traceless and singular 
exact solution 

X=M, with trM = trM^ = detM = , (66) 
where M is any 4x4 constant matrix with trM = trM^ = detM = and arbitrary trM^ & rank (M). 

Corollary 6. 

Any traceless solution to (ITSl ) & (fT9l ) is a singular constant matrix (l66l ) with / = 0, ^ = 0, ^ = and arbitrary h. 

Looking for constant solutions to the system ([TSl l & ([T9l ). the solution (l66l ) was derived analytically in 
Ref. [3], but it has not been acknowledged in the literature that it is the unique traceless solution to ([TSl l & (fT9l) . 

7 Conclusion 

Using ADecM we have shown, relying on the particular case ([TS] ). that solutions x{p) to equation Q are 
necessarily polynomials in a constant matrix with scalar coefficients, which are functions of p (equation 
Restricting ourselves to the case of a 4x4 nonlinear MDE along with a scalar one, which constitute the field 
equations of a 5DGR extended by a quadratic GB term in the SCS case (equations ([TSl l & ([T9l)). we have been 
able to construct the general unique power series solutions (equations (|49l ) & (|54l)). These power series terminate 
in closed-form singular (det % = 0), nontraceless (l52l) and traceless (l66l ) solutions, which all satisfy the property 
tr;^^ — (tr^)^ = 0- No nonsingular (det;^ / 0) or further singular solutions have emerged from the unique power 
series solutions. This constitutes a proof of uniqueness of the re-derived solutions (l52l) & (l66l) . 

In summary, we have developed in this paper an efficient ADecM to solve the constrained nonlinear MDE's 
of type (dH). The search for SCS solutions in JDGR {d > 5) leads to a system of {d - \ ) x {d - \) MDE along 
with a scalar one, solutions of which may be laid down by the ADecM developed here. 

Acknowledgments 

The author wishes to thank Gerard Clement (LAPTH) for useful comments on this manuscript. 

Appendix A: Proof of Corollary 3. 

In equation ( |49l ) we set Mi = 2A. The constraints on A become: trMi = 2 & trM^ = 4. Now, ^ being provided 
by ( |49l ). the general expressions up to the order 7 for G and k read 

G = [trM? - {tvMif]x^ + [2tr (M1M3) - 2trMitrM3]/ 

+ [2tr (M1M5) + ti-Mf - {tvM^f - 2ti-MitrM5]/ + ff{x^) , 
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k = [(trMi)'* - 6{trMi)hrMj + 3{trM^f + StrMitrM^ - 6trM/]//24 

+ [2(trMi)^trM3 - GtrMitrM-itiMj - 6(trMi)2tr (M1M3) - 12tr (Mj^Ms) 
+ 6trMftr(MiM3) +4trM3trMj' + 12trMitr(MfM3)]//12 + ^(x^) , 

Substituting the above expressions and values into dTSl ) & ( fT9l ) and using the Cayley-Hamilton equation for Mi, 
the terms of order 4 in these equations lead to the matrix & scalar equations, respectively 

47(2M? -Mi^) +2(trM3)Mi - 8M3 - 16trM3 +4tr (M1M3) =0, (A.l) 
37detMi +4trM3 - 2tr (M1M3) =0. (A.2) 

In general the terms of order l + l (l odd) in (fTSl ) & ([T9l ) lead to matrix & scalar equations, respectively, depend- 
ing linearly on trM/ & tr (MiM/) as in (lA.ll ) & (IA.2I) . Furthermore, the matrix equation depends linearly on M/. 
Such a system can always be solved for (trM/, tr(MiM/), M/) as we are going to do for (lA.ll) & (IA.2I) . Since 
the resolution of the system involves tracing the matrix equation, consistency of the obtained solution (trM/, 
tr (MiM/), M/) has to be checked for each step. 

Solving the system consisted of (IA.2I ) and the trace of (lA.ll) . we obtain 

2 y 

trM3 = -7detMi + ^(8-trMf) (A.3) 

tr (M1M3) = ^ydetMi + ^(8 - trM^) . (A.4) 
6 9 

Substituting these last two equations in (lA.lb . we obtain 

M3 = |(2M2 -Mi^) + i (^^detMi + |(8 - trM3)^Mi 

+ XdetMi-|(8-trM3). (A.5) 

Tracing (IA.5I ) reduces to (IA.3I) and tracing (IA.5l) xMi reduces to (IA.4I ). Hence, the solution (trM3, tr(MiM3), 
M3) is consistent. 

We applied the same procedure for the order 6 and obtained the following solution (trMj, tr (MiM5),M5) 

trMs = ^ [1024 + 1698 detMi + 171(detMi)^] 

- (1287+967detMi) + ^^'^j^^'^' , (A.6) 

72[-512 + 13164 detMi+2223(detMi)2] 



tr(MiM5) 

y('-128y+ 1929ydetM, ItrM? v^itrM?)'^ 

(A.7) 



1944 

7(-1287+ 19297detMi)trM3 Y^itrMfy 
1944 243 
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+ 1-9-+ 48 ^^r' 

254/2 54772detMi 3l72(detMi)2 73T^trMf 
243 ^ 1296 ^ 864 486 
1972detMitrM^ ^ \9f-{\xM\Y\^ ^ ^-fidi^mxf ^if-dL^VMi 



648 7776 y 48 72 

1372trM^ 5272 72detMitrMj' \^f{V\:MXf 
^ 27 27~"^ 9 432 ' 

The solution (trMs, tr(MiM5), M5) failed to be consistent. In fact, subtracting (IA.6I ) from the trace of (IA.8I ) and 
subtracting (IA.7I ) from the trace of (IA.8l) xMi lead to the following two consti^aints on the two free parameters 
(trM^, detMi) 

457(detMi)2-137(trMi3-8)2 + 67detMi(8trM^-37) =0, (A.9) 
907(detMi )2 - 87(trM^ - 8)^ + 37detMi (89trM^ - 496) = . (A. 10) 

Substituting the series expansion of % (equation (l49l )) into its Cayley-Hamilton equation ([T5] ) and replac- 
ing the matrices M3 , M5 and their traces by the expressions derived above without, however, using the con- 
straints (IA.9I ) & (lA.lOl ). we obtain upon rearranging the terms the following "matrix" series expansion 

{[(167 + 67detMi+9(trM2)^)trM^-647-487detMi -7(trMj')2]Mi (A.ll) 
+ 967detMi - 1287+(327- 127detMi)trMj' - 27(trMi')2}xV36 + ^(x^) = 0. 

Notice that the matrix coefficients of & in (lA.l II ) are identically zero. Since (lA. Ill ) is an identity, all the 
matrix coefficients have to be zero. Hence, tracing the matrix coefficient of x^ in (lA. Ill ), we obtain the following 
constraint 

[187detMi + 57(trMi' - 8)] (trMj^ - 8) = . (A. 12) 

Solving the system of the three constraints ( |A.9b . (I A. 101 ) & (IA.12I ) leads to 

(detMi = & trMf = 8) ^ (detA = & trA^ = 1) . 



Appendix B: Proof of Corollary 6. 

Given / = 0, trMo = & trM^ = 0, we have shown in subsection I6.2l that trMp & Mq remain arbitrary and that 
detMo = & Ml = 0. To show that the other matrices Mj, j > 2, are also zero we conduct a proof by induction. 
Knowing that Mi = we assume that, for 7 > 2, Mi = M2 = • • • = M^_i = and will show that Mj = 0. X' 
provided by (l54l ). and its m"' power become 

X = Mo + MjP^ + Mj+i pj+^ + ■■■, 
With Ml = M2 = • • • = M^_i = for j > 2, equations ([55]) become 

00 00 00 

£ C„p" = 0,([l9l)^ £c„p" = 0,(115]) ^ ££„p" = 0. 

n=j- 1 n=j n=j 
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Since for / = we have tr Mj = 0, the scalar equation cj = leads immediately to tr (MoMj) = 0. Using this 
last equation in the matrix equation Cj-i = we reduce it to 

{4 + 6YM^)Mj = 0. (B.l) 

Multiplying both sides of (IB.ll ) by (1/4) — S/Mq /8, which is the inverse matrix of 4 + 6yMq, leads to Mj = 0. 
With this value of Mj, the matrix equation Ej = vanishes identically. 
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